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Abstract. 

Himalayan mountain belt was created by the Indo—Eurasian collision. The length of 
this seismically active mountain belt is almost 2500 km. Southern partition of this zone is 
matked by the main frontal thrust. We computed the isobase, relative relief, incision, vertical 
dissection and drainage density maps of the study site. Drainage density is inversely 
proportional to the rest of computed parameters. Drainage density is less where the 
incision/deep cutting is high. Pezu and north western parts of the study site are observed 
vulnerable to tectonic activity and high risks. These sites are the adjoining parts of active 
faults. The rates of morphological changes and the combination of stripping processes are 
determined through Surface dynamics maps (SDM). Although it is a handy tool used to 
evaluate erosion conditions of sags or sub basins. Remote sensing and GIS techniques 
proved efficient for appraisal of thrust in Marwat Kishore and Bhittani ranges. 
Keywords: Thrusting, Relative Relief, Isobase coding, Drainage density. 
Introduction 

The Himalayan mountain belt was created by the Indo—Eurasian collision. The 
length of this seismically active mountain belt is almost 2500 km. Southern partition of this 
zone is marked by the Main Frontal Thrust (MFT) [1,2]. The growth and development of 
MFT is necessary to unwind the process of deformation of mountain belts and preparation 
of major jeopardy of future related to active tectonics of these belts. The data gained 
through geodesy suggests that the Indian plate move at the rate of 35mm/year in 
northwestern boundary while 38mm /year in northeastern direction [3]. According to Lave 
and Avouac most of the Himalayan deformations were recorded in MFT [4]. The 
deformation was deliberated 9 + 3 mm/year in the northwestern, 21 + 1.5 mm/year in 
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central and 23.4 + 6.2 mm/year in eastern Himalayas [4,5,6,7,8]. The MFT is a youngest 
frontal thrust and a narrow region located in the central and eastern Himalayas which has a 
high degree of cross-section; however, a low degree of cross-sectional taper is located in the 
western Himalayas which is 100 km wide [9,10]. Moreover, the Indian plate rotates in an 
anticlockwise direction shifting the trend from NW-SE of central and eastern Himalayans to 
NE-SW located in Pakistan [11]. Himalayan arc is mapped in Figure 1. 
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Figure 1. Himalayan arc observed through satellite imagery 
Soutce:https://www.mdpi.com/2072-4292/12/20/3362/htm. 

Several historic earthquakes which had a magnitude higher than seven were 
associated with the MFT. These earthquakes had the ability to rupture the surface of land. 
The Himalayans located in the north of Pakistan elongates from Jhelum fault in east to 
western Ophiolite Belt. The main frontal thrust of Western Himalia is classified as salt Indus 
ranges and trans Indus ranges. 

According to geodetic data, southward translation of western salt rang is 10mm /year 
and central salt range is 5mm/year. However, Kalabaghfault zone displaced the Trans Indus 
Range (TIR) at the rate of 5.3 mm/yeatr. Blisniuk et al. and Abir et al. [12] suggested that the 
basal detachment is formed along the MFT by the salt range in the west of Himalayas. 
Chronostratigraphic studies proposed that western Himalayans are youngest collisional zone 
which is deforming actively. 

The trans indus ranges are further sub-divided into the Surghar, Manzai, Marwat— 
Khisor, and Bhittani Ranges (Figure 2). These ranges portray condensational deformation 
and show lateral structural variations along the salient faults. Moreover, transpersonal 
deformation occurred along strike-slip faults [10]. The Kalabagh and Tank faults are 
separated by Marwat—Khisor salient in the TIR. 
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The main objective of this study is to determine the deformation pattern, 
displacement rate and landscape development relative to active tectonics of western 
Himalayans. 
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Figure 2. Fault setting in the study stie. 
MATERIALS AND METHODS 
Study Area 


The atea under study covers about 7850 square kilometers of Bannu, Kalabagh, 
Lakki Marwat, Pezu, Tank and D.I.Khan (Figure 3).The study area is situated about 200 
Kilometers away from Peshawar towards south-southwest. 
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Figure 3. Ranges present in the study site 


Data Acquisition and Processing 

Data was obtained through geological fieldwork for geological mapping of trans 
Indus ranges and geological cross-sections were established to comprehend geometry and 
structure of Trans Indus range. Two types of data were used for active tectonic mapping 
which include vector data and Shutter Radar Topographic Mission (SRTM) Digital Elevation 
Model (DEM) [13]. 

Digital Elevation Model (DEM) 

DEMs ate helpful in determining geographic information based on digital relief 
maps. However, DSM (Digital Surface Model) was used in landscape mapping and city 
modeling. 

Longitudinal River Profile analysis 

A flow direction of any topographic region is determined through D8 algorithm. In 
this method, flow from one pixel to its neighboring pixels is tracked. D8 flow grid is 
generated from a depression less DEM. This method is used to determine flow from flat to 
reducing valleys. Incision determine the nature of geology, soft rocks determine erosion 
while hard rocks lead to vertical cutting/incision. Relative relief describes the variation in 
slop, high relative relief leads to gentle slopes while low relative relief describes the steep 
slopes. Drainage density was computed using hydrological tools on digital elevation model, 
where large drainage density leads to rough topography. 

Concavity and steepness 

Stream power model is used to estimate the drainage area and channel gradient 

(Whipple and Kirby, 2001; Kerby and Howard 1983): 
dz! dt =U(x,t)— KA"S" fq. 1 

Where dz/dt is the rate of change in variation of channel’s elevation, U determine 
the rate of uplift in rock, A describe the area of drainage, S is the slope of each channel, K is 
the coefficient of erosion and m and n describes the hydrology and hydraulic geometry 


respectively. The value of dz/dt is zero in steady state conditions. To compute the slope of 
channel, following expression is used, 


5 _ (U / Ky" Aen) Eq. > 


Equation 2 describes a relationship between the area of drainage line and the channel 
gradient which is commonly observed in natural conditions. 


y= kA“ Eq. 3 


Here the factor (U/IX)1/n determine the steepness, ks and m/n is the concavity. We 
simply applied log on both sides of equation 3, 


log S =—@log A+logk, 44.4 


Where 0 determine the concavity and K, is the steepness which may be computed by 
regression in log plot of slope and drainage area. 
RESULTS AND DISCUSSIONS 
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Hack gradient index and geomorphic indices were analyzed in the study site in order 
to interpret relative uplifting, regional topography and crustal slanting within the region. The 
concavity (0) and steepness (ks) indices were calculated using stream power law by 
establishing stream longitudinal profile analysis on 116 small and large channels which were 
extracted through the SRTM DEM. We applied regression models to active tectonic signals 
from streams in order to determine concavity indices and steepness. The steepness index 
values fluctuate between high or low concavity. Spatially inconstant lithologies, or spatially 


variable relative bedrock uplift rates causes such a fluctuating zone. 
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Figure 4. Drainage Density Map at 4Km window. 
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Figure 5. Vertical dissection window at 4Km window. 
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Figure 6. Incision Map at 4Km window. 
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Figure 7. Relative Relief Map at 4Km window. 
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Figure 8. Isobase Coding Map at 4Km window. 
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Figures 4-8 determine the isobase, relative relief, incision, vertical dissection and 
drainage density maps of the study site. Drainage density is inversely proportional to the rest 
of computed parameters. Drainage density is less where the incision/deep cutting is high. 
Pezu and north western parts of the study site are observed vulnerable to tectonic activity 
and high risks. These sites are the adjoining parts of active faults. 

The rates of morphological changes and the combination of stripping processes are 
determined through Surface dynamics maps (SDM). Although it is a handy tool used to 
evaluate erosion conditions of sags or sub basins. Differential erosion and neotectonics 
activity were found behind the uplift of study site due to abundant precipitation in summer 
as well as in winter and tectonic domination over climate and erosion. 
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